6-28  quad8ʹӦ(Simpson)֡
 
ȱдMļ¡
function f=fun2(x)
f=x.^2./exp(-x);
end
Enterȷϡ
>> quad8('fun2',0,3)
ans =
   98.4277
>> a=quad8('fun2',1,3,[1e-10,1e-11])
a =
  97.7094
>> vpa(a,10)
ans =
97.70940279
>> quad8('fun2',1,3,1e-6,1)
      18     1.0000000000    1.00000000e+000    97.7094029629
ans =
   97.7094
>>
